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Abstract 



We present an analytic study of the mode-coupling phenomena for a scalar 
field propagating on a rotating Kerr background. Physically, this phenomena 



> 

(N 
(^ . 

^^ , is caused by the dragging of reference frames, due to the black-hole (or star's 

^^ \ rotation. We find that different modes become mixed during the evolution 



and the asymptotic late-time tails are dominated by a mode which, in general, 
has an angular distribution different from the original one. We show that a 
rotating Kerr black hole becomes 'bald' slower than a spherically-symmetric 
Schwarzschild black hole. 

I. INTRODUCTION 

The radiative tail of gravitational collapse decays with time leaving behind a Kerr- 
Newman black hole characterized solely by the black-hole mass, charge and angular- 
momentum. This is the essence of the no-hair conjecture, introduced by Wheeler in the 
early 1970s |l|. The relaxation process of neutral fields was first analyzed by Price [§ for 
a nearly spherical collapse. Price demonstrated that the fields decay asymptotically as an 
inverse power of time. Physically, these inverse power-law tails are associated with the 
backscattering of waves off the effective curvature potential at asymptotically far regions 



The existence of these inverse power-law tails, which characterize the asymptotic late- 
time evolution of test fields on a Schwarzschild background was later verified by other authors 
|Q-|^ using independent analyses. The late-time evolution of natural test fields on different 
spherically symmetric spacetimes was studied in Refs. 0,0,0 • In addition, the late-time tails 



of gravitational collapse were demonstrated by full numerical computations |10|,|TI|. The 
physical mechanism by which a charged black hole, which is formed during a gravitational 
collapse of a charged matter, dynamically sheds its charged hair was first studied analytically 
in [p]2| , p!3| , and a full non-linear numerical analysis was performed in |l^ . 



However, all the above mentioned analyses were restricted to spherically symmetric back- 
grounds (in particular, to the Schwarschild and Reissner- Nordstrom black holes). On the 
other hand, astrophysical objects are usually rotating. Hence, the physical process of stellar 
core collapse to form a black hole is expected to be non-spherical in nature. The corre- 
sponding problem of wave dynamics outside a realistic, rotating black hole is much more 
complicated due to the lack of spherical symmetry. Originally, this problem was addressed 
numerically by Krivan et. al. |T^|T^, followed by Ori |T^ and Barack [|T^ who provided 
preliminary analytical results for the evolution of a test scalar field on a Kerr background. 

Evidently, the most interesting situation from a physical point of view is the propagation 
of gravitational waves on a rotating Kerr background. The evolution of higher-spin fields 
(and in particular gravitational perturbations) outside a realistic, rotating black hole was 
studied analytically only recently |19[ . This was done by analyzing the asymptotic late-time 



solutions of Teukolsky's master equation p0|j21[| , which governs the evolution of gravitational, 
electromagnetic, neutrino and scalar perturbations fields on Kerr spacetimes. 

The main result presented in Ref. [|1^] is the existence of inverse power-law tails at the 
asymptotic regions of timelike infinity z+, null infinity scri^ and at the black- hole outer 
horizon iJ+ (where the power-law behaviour is multiplied by an oscillatory term, caused 
by the dragging of reference frames at the event horizon). This late-time behaviour is 
similar to the late-time behaviour of massless fields on a spherically-syTaTaetnc Schwarzschild 



background. 

However, there is one important feature (besides the oscillatory behaviour along the 
black-hole horizon) which is unique to rotating collapse: Active coupling of different modes. 
The phenomena of mode-coupling has been observed in numerical solutions of Teukolsky's 
equation |TB|,|TB[ (a similar phenomena is known in the case of rotating stars |2^). Physically, 



this phenomena is caused by the dragging of reference frames, due to the rotation of the black 
hole (or the star). The issue of mode-coupling was discussed in Ref. [|l9l, where preliminary 
considerations and results were presented. In this paper we give a complete mathematical 
analysis of this important phenomena for the case of a scalar field propagating on a realistic 
Kerr background, and we focus on its physical implications. We consider all the various 
possible cases of mode-coupling, which result in different types of asymptotic late-time 
evolutions (different damping exponents). 

We find it most convenient to use the basis of the (frequency-dependent) spheroidal wave 
functions and carry out the analysis in the frequency domain; these functions allow us to 
separate the ^-dependence from the radial dependence, and, more important, allow us to 
separate the different modes (different values of /) so each mode evolves independently. On 
the other hand, our main objective is to study the temporal dependence of the field's late- 
time behaviour; Therefore, at the end of the analysis we should integrate over all frequencies 
in the frequency plane to obtain the explicit temporal dependence of the field. Thus, the final 
expression for the field is most naturally expressed in terms of the u- independent spherical 
harmonics. 

The plan of the paper is as follows. In Sec. |I| we give a short description of the physical 
system and summarize the main analytical results presented in Ref. |19|. In Sec. |T| we 
discuss the effects of rotation and the mathematical tools needed for the physical analysis 
are derived. In Sec. |rV| we study the active coupling of different modes during a rotating 
gravitational collapse, with pure initial data. In Sec. |V| we consider the coupling of different 
modes, with generic initial data. In the appendix we study the coupling of the real and the 
imaginary parts of the radial field. We conclude in Sec. ^ with a summary of our analytical 



results. 



II. REVIEW OF RECENT ANALYTICAL RESULTS 



We consider the evolution of a massless scalar field outside a rotating star or a black hole. 
The external gravitational field of a rotating object of mass M and angular-momentum per 
unit-mass a is given by the Kerr metric, which in Boyer-Lindquist coordinates takes the 
form 



ds' 



2Mr\ ., 
1 T^\dt' ^ 



sin^ e\r'^ \a^ \ 



2 , ^4Marsin2^^|^^^^ ^ ^^^ ^..2 
2Ma^rsin^e' 



dtdip - —dr'^ - T.de'^ 



d^' 



(1) 



where T, = r'^ + a? cos^ 9 and A = r^ — 2Mr + 0? (we use gravitational units in which 
G = c=l). 



In Boyer-Lindquist coordinates the wave equation for a massless scalar field reads ||23 



{P + a^f 
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Resolving the field in the form 



E ^' 



(3) 



where m is the azimuthal number, one obtains a wave-equation for each value of m: 



d^^ (9* d'^^ A 1 d f d^\ 

B,(r,0)—^B,(r.S)-- — ^B,ir,0)^-^^^^^ — -[.^0-y, , (4) 

where the tortoise radial coordinate y is defined by dy = ^ ~^ dr. The coefficients Bi{r,9) 
are given by 

Aa^ sin^ 9 



Bi{r,9) = l 



[r^ + a 



2^2 



(5) 



B.(r) 



AiMmar 



(6) 



and 



Ssfr, 



2(r-M)r(r^ + a^) 
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771 



A sin^ e 
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(7) 



The time-evolution of a wave-field described by Eq. (P is given by 



^(z,t) = 27r 



5i(z') G'(^, ^'; t)^i(/, 0) + Gi(z, /; t)^{z', 0) 



+ 



52(2')G'(^>-2';^)*(^',0)Uin^Wrfy' 



(8) 



for t > 0, where z stands for {y, 6), and G(z, z'; t) is the (retarded) Green's function, which 
is defined by 



B,(r.»)^ + ft(r.<^)|-|, + B3(r. 



A I d ( . d^ 

(r^ + a^y sine m v^^ m 



G{z,z';t) 



S{t)S{y-y 



.5(6-9') 



(9) 



2tt sin 9 

The Green's function (which is associated with the existence of a branch cut in the complex 
frequency plane 0, usually placed along the negative imaginary u— axis) was calculated in 
Ref. ^ [see Eq. (31) of [|l|]]: 

.1)1+1221 ^-i^ r(/ + 1 - 27ujM) 

(10) 



^ ^^'^'^^ ,^|7rA2[(2/ + l)!]2yo T{-l-2iu;M) 

^i{y, uj)^i{y', uj)Si{9, auj)Si{9', auj)uj^^+^e~''^'duj 



where the function '^i{y',uj) is given by (see Ref. [T^ for additional details): 



^1 = Ar'+'e*'^^M(/ + 1 - 2iLuM, 21 + 2, ~2iujr) 



(11) 



and A is a normalization constant. This expression for the Green's function is obtained 
using the large-r (or equivalently, the law-cj approximation ||19|). The functions Sp{9,auj) 
are the spheroidal wave functions which are solutions to the angular equation ||2^j2T[ 



1 ^ ^ ■ n^\ , 2 2 2n 

sin9-—- \ + a 00 cos 9 



m 



+ AT 



^r = . 



(12) 



_sin9 89 y"^ " 39 ) ' " ^™ " sin^^ 
For the ao; = case, the eigenfunctions SY^{9,auj) reduce to the well-known spheri- 
cal harmonics Y{^{9,(f)) = S'™(0)e*™'^, and the separation constants AJ^lauj) are simply 



III. ROTATION EFFECTS 

The rotational dragging of reference frames, due to the rotation of the black hole (or 
the star) produces an active coupling (during the gravitational collapse) between different 
modes. Thus, even if the initial data is made of a pure mode (i.e., it is characterized by a 
particular spherical harmonic 1^™) other modes would be generated dynamically during the 
evolution. Mathematically, the coupling of different modes is attributed to the 6'-dependence 
of the spheroidal wave functions Sl^{6, auj) and of the coefficient -Bi(r, 9). 

It is well known that the late-time behaviour of massless perturbations fields is deter- 
mined by the backscattering from asymptotically far regions P,@]. Thus, the late-time 
behaviour is dominated by the /oty-frequencies contribution to the Green's function, for only 
low frequencies will be backscattered by the small effective curvature potential (for r ^ M). 
Thus, as long as the observer is situated far away from the black hole and the initial data 
has a considerable support only far away from the black hole, a small-uo approximation is 
sufficient in order to study the asymptotic late-time behaviour of the field 0. 

If the angular equation (p!2D is written in the form 

(L° + L^)S]f' = -A^Sl^ , (13) 

with 

and 

L\e, auj) = {aujf cos'^e , (15) 

it is immediately amenable to a perturbation treatment for small au. The spherical functions 
Yf^ are used as a representation. They satisfy L^Yi = —A\ Yi with A\ = /(/ + 1) (we 
suppress the index m on Ai and Yj). For small au standard perturbation theory yields (see. 



for example, p5|) 



\L 



1 |2 



and 



s, = Y, + y: .,„, " .,„, yt + 



A, = A'°'+L; + E-;TiFV + -- (16) 
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f^ 4(0) _ 4(0)- 



k^lnjLl[Ai - Af^ }{Ai -An ) k^l[Ai - Af, } 



+ ■■■ (17) 



where 



Ll = {Olm\L^\Okm) = f Yi"'*L^Y^'dQ . (18) 

Equation ([17|) implies that the black-hole rotation mixes various spherical harmonics func- 
tions. 

The spherical harmonics are related to the rotation matrix elements of quantum mechan- 
ics p6[. Hence, standard formulae are available for integrating the product of three such 
functions. The one we need in order to evaluate the integral Eq. ([T8|) is 



2. 1., V 1. , 2 f2k + l 



1/2 



(0/m| cos^ e\Okm) = -6^ + - , (A;2mO|/m)(A;200|/0) , (19) 

where {jiJ2fnifn2\JM) is a Clebsch-Gordan coefficient. The well-known properties of the 



Clebsch-Gordan coefficients 27 : 



1 • ( ji J2"^i"^2 1 JM) vanishes unless ji , J2 and J satisfy the triangle selection rule | ji — J2 1 < 

J <jl+J2- 

2. (J1J2OOIJO) = if ji + J2 + J is an odd integer (this is a parity conservation law). 

imply that Lji^ 7^ for / = k,k ±2 (all other matrix elements vanish). In general, if we 
choose the spherical harmonics as a representation, so that 

00 
Sii9,au;) = ^ Gfc(aa;)l'-'^ln(e) , (20) 

k=\'m\ 

then, to leading order in (atu)^, the coefficients Cik{auj) take the form (see, for example, p5| ) 



k-2 



{On + 2m\ cos^ 6\0nm) 



Cik - n — 7(^) — jm — ' (21) 



/I vu; Ay 

n=l,l+2,... ^l — ^n+2 

a k > I + 2 and k — I is even, and 

t^ {Onm\ cos^ 9\0n + 2m) 
^ik = 11 -joj -(0) . (22) 

n=fc,fc+2,... ^i ~ ^n 

a k < I — 2 and / — fc is even. In addition, C/^ = if |/ — /i;| is odd, and Cu = 1. 

The time-evolution of a wave-field described by Eq. (Q) is given by Eq. @. The 
coefficient -Bi(r, 6*) appearing in Eqs. (^ and (|) depend explicitly on the angular variable 9 
through the rotation of the black hole (no such dependence exist in the non-rotating a = 
case). To elucidate the coupling between different modes we should evaluate the integral 
[see Eqs. (H) and (|)] 



2.,^, V 2, 2 (2k + l 



1/2 



{Olm\ sm^e\Okm) = -Ski - - —, (A;2m0|/m)(fc200|/0) . (23) 

Clearly, this integral vanishes unless I = k,k ± 2. This, togather with the expressions for 
the coefficients Cik imply that if a mode YiT' presents in the initial data, then all modes Y^ 
with |/* — /| even (and I > \ni\) will be present at late times. 

As long as we are interested in the asymptotic late-time behaviour (the u ^ limit) of 
the solution, we may use the approximation 

T{-l-2ttuM) V ; ' y J 

which is valid for tu ^ 0. Thus, we obtain [see Eq. (|ToD] 



m 
—ioo 



7rA2[(2/ + l)!!]2 
§i(y, uj)^i{y', uj)Si(e, auj)Si(e', acu)cj2'+2e-^-*rfu; . (25) 

A note is needed here on the separation of the r and 9 variables in the Kerr background. 
These variables are separated by Teukolsky's equation only in the frequency domain. This is 
caused by the fact that the spheroidal harmonics which are used for the separation depend 
explicitly on the temporal frequency uj. On the other hand, our main objective is to obtain 

8 



expressions for the late-time behaviour of the field in terms of the temporal variables (t, u, 
and f ); at the end of the analysis we should recover the explicit temporal dependence of the 
field by integrating over all frequencies in the frequency plane (this would lead to the final 
answer which is, of coarse, cu- independent). 

For a generic initial data in the time domain, one is able to write these 

data as a sum over l,m, and integration over u. Namely, ip{r,6,ip,t) = 

E E e'"''^ JZoRrir,t)Sl^{e,auj)e-''^^du} with t = 0. Then, each mode [characterized 

m=—oo l=—m 

by a particular spheroidal wave function Sp{6,auj)] evolves independently. Nevertheless, 
Eq. (pOD implies that each spheroidal wave function can be expressed as a sum over the (a;- 
independent) spherical harmonics (were the cj-dependence of the spheroidal wave function 
is now moved to the expansion coefficients). That is, the initial data can be expressed also 

oo m 

as -^{r, e,!f,t)= E E Rr{r, t)Y{^{e)e''^'^ with t = 0. 

m=—oo l=—m 

It is convenient to use the basis of the (tu-dependent) spheroidal wave functions and 
carry out the analysis in the frequency domain; these functions allow us to separate the 
^-dependence from the radial dependence, and, more important, allow us to separate the 
different modes (different values of /) so each mode evolves independently, see Eq. (12) of 
|IP|. However, since our goal is to calculate the temporal dependence of the field [in terms 



of the time t (or u, v)], the final expression for the field should not involve the u— dependent 
spheroidal harmonics. Therefore, the final expression for the field is most naturally expressed 
in terms of the uj- independent spherical harmonics. 

We shall consider two kinds of initial data: A generic initial data, which can be decom- 
posite in terms of the basis of the spheroidal wave functions. Then each spheroidal wave 
function evolves independently. In view of Eq. (0), this corresponds to the assumption that 
the (generic) initial pulse consists of all allowed modes (all spherical harmonics functions 
with / > \m\). 

On the other hand, in order to compare our analytical results with the numerical results 
of Krivan et. al. [|r5[, one should consider an initial data which is made of a single mode. 



This corresponds to the assumption that the initial angular distribution is characterized by a 



particular spherical harmonic function Y^T". We begin with this case, which can be compared 



with the numerical results of 115 



IV. PURE INITIAL DATA 

A. Asymptotic behaviour at timelike infinity 

We consider first the behaviour of the scalar field at the asymptotic region of timelike 
infinity i+. As was explained, the late-time behaviour of the field should follow from the 
/oiy-frequency contribution to the Green's function. Actually, it is easy to verify that the 
effective contribution to the integral in Eq. (|25|) should come from |c(j|=0(j). Thus, in order 
to obtain the asymptotic behaviour of the field at timelike infinity (where y, y' <ti t), we may 
use the \uj\r ^ 1 limit of \E'i(r, cj). Using Eq. 13.5.5 of ^^ one finds 

^i(r,cu) ~ylr'+^ . (26) 

Substituting this in Eq. (^) we obtain 

qC( h) ^ y- j^ (vv'Y+^ 



r—ioo 

/ Siie,auj)Si{e',aiu)uj^^-^^e-''^'duj . (27) 

Jo 



Using the representation for the spheroidal wave functions Si, Eq. (pOD, togather with 
the fact that the integral Eq. (|23| ) vanishes unless I = k,k ±2, we find that the asymptotic 
late-time behaviour of the / mode (where |/* — /| is even and / > \m\) is dominated by the 
following effective Green's function: 

^ 9'iM 

/•—'too 

CkiCki*^M0)Yi:_,i9') / (a^)'*+^-2'=-'^a;2fc+2g-.c.t^^ ^ (28) 

'J u 

where p = if /* — \m\ is even and p = 1 otherwise, and g = 2 if /* > |r?7,| + 2 and g = 
otherwise. Here, L = I* — q for / > /* modes, and L = I ioT I < I* — 2 modes. Performing 
the integration, we obtain 
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k=\m\+p 



TX 



[(2fc + l)!!]2 



Thus, the late-time behaviour of a scalar field at the asymptotic region of timelike infinity 
i+ is dominated by the lowest allowed mode, i.e., by the I = \m\ mode if I* — \m\ is even, 
and by the I = \m\ + 1 mode if Z* — \m\ is odd. The corresponding damping exponent is 
-(/* + \m\ +P + 1) if /* > \m\ +2, and -(2/* + 3) HI* = \m\,\m\ + 1. 

These analytical results should be compared with the nwmenca/ results of Krivan et. al. 



||15|| . The comparison can be made for the specific initial data considered in ||T5[, and is 
presented in Table ^ - we find a perfect agreement between our analytical results and the 
numerical results of Krivan et. al. [|1^ . 

B. Asymptotic behaviour at future null infinity 

We further consider the behaviour of the scalar field at the asymptotic region of future 
null infinity scri^. It is easy to verify that for this case the effective frequencies contributing 
to the integral in Eq. (125|) are of order O(-). Thus, for y — y' <^t <^2y — y' one may use 



the \uj\y' <^ 1 asymptotic limit for \l'i(y', cj) and the \uj\y ^ 1 {Imuj < 0) asymptotic limit 
of \^i(y,u;). Thus, 

^i(y',a;)^ V\ (30) 

and 

^i{y,u:) ~ Ae-y{2l + 1)!^ ^,, ^ ' -— ^ , (31) 



where we have used Eqs. 13.5.5 and 13.5.1 of |^^, respectively. Substituting this in Eq. 
(P5|) we obtain (using the Muo -^ limit) 

^c. ,.,. _ f M2-'(2/ + l)!(-l)('+^)/^ 
^''^"il vr/![(2/ + l)!!p y 



-lOO 







Si{d,auj)Si{d',au)J+^e-''''^'-y^duj . (32) 
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Using the representation Eq. (pO|) for the spheroidal wave functions 5*/, togather with the 
fact that the integral Eq. (|23| ) vanishes unless / = k,k±2, we find that the behaviour of the 
I mode (where |/* — /| is even and / > |m|) at the asymptotic region of null infinity scri+ is 
dominated by the following effective Green's function: 

' ^ ' ' ^ " .>_,, nki[i2k+m^ y 

r—ioo 

CkiYi{e)Y*{9') / (acj)'-'=cu'=+^e-*^(*-^)rfw , (33) 

Jo 

where 52 = for the / = /* mode and ^2 = 2 for I > I* + 2 modes (gi = 2 if /* > \m\ + 2 and 
gi = otherwise), and 

.c,_^.^_M2-'(2/ + l)!(-l)('+^)/^.+l 



Gf{z,z'-t) 



7r/![(2/ + l)!!]2 

r—ioo 

Ca*-2Yi{e)Y;_,i9') / iacuf'^-'J+'e-^'^^'-y^diu , (34) 

Jo 



for I < I* — 2 modes. Performing the integration in Eqs. (p3D and (|3^), one finds 

^*+^^ M2-^'(-l)('+^+2)/2(2A: + l)!(/ + l)! ,,^, 

' ^ ' ' ^ " .i,. .ki[i2k+m]^ y 

CMYi{e)Y:{ey-'u~'^'+''> , (35) 

for / > /* modes, and 

_ M2-H-l)('*+')/^(r - 1)1(2/ + 1)! ,^, 
^'^'^''^'*^" 7r/![(2/ + l)!!p y 

Cu,.2Yi{9)Yi:_^{9')a''-'-'u-'* , (36) 

for I < r — 2 modes. Thus, we find that the late-time behaviour of a scalar field at the 
asymptotic region of null infinity scrij^ is dominated by the I < I* — 2 modes provided that 
/* > \m\ + 2 (where \l* — l\ is even and / > Iml). These modes decay as an inverse power of 
the retarded time u, with a damping exponent equals to — /*. If /* = \m\, \m\ + 1 then the 
late-time evolution will be dominated by the I = I* mode, which has a damping exponent 
equals to -(/* + 2). 
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C. Asymptotic behaviour at the black- hole outer horizon 

Finally, we consider the behaviour of the field at the black- hole outer- horizon iJ+. The 
asymptotic solution at the horizon {y — >• — oo) is given by |^ (see Ref. [T^ for additional 
details) 

§i(y, to) = c(cu)e-'("'-™"'+)^ , (37) 



where w+ = a/(2Mr+). In addition, we use Eq. ( pOf ) for "^1(1/', lu). 

Using the representation Eq. ( PP] ) for the spheroidal wave functions Si, togather with 

the fact that the integral Eq. ( PB| ) vanishes unless I = k,k ±2, we find that the asymptotic 

behaviour [where C{uj) = const. + 0{Muj)] of the / mode (where \l* — 1\ is even and / > \m\) 

at the black-hole outer horizon H^ is dominated by the following effective Green's function: 

rC(,,'.t)- f 2M(-l)('-+^+^-^)/^(r + / + 2-g)! 

where q,p and L are defined as before, and Tk are constants. Hence, we find that the late- 
time behaviour of a scalar field at the black-hole outer horizon is dominated by the lowest 
allowed mode, i.e., by the / = \m\ mode if /* — |m| is even, and by the / = \m\ + 1 mode if 
/* — \m\ is odd. 

V. GENERIC INITIAL DATA 

In this section we consider the generic case. That is, we assume that the initial pulse 
consists of all the allowed (/ > \m\) modes (see also the recent analysis of Barack and Ori 
p8|). The analysis here is very similar to the one presented in fV[ We consider first the 
behaviour of the scalar field at the asymptotic region of timelike infinity i+. Using Eq. (p7D, 
togather with the representation Eq. ( pOf ) for the spheroidal wave functions, we find that the 
asymptotic late-time behaviour of the / mode (where / > |r7i|) is dominated by the following 
effective Green's function: 
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r — ioo 

CMW)Y:{e') / (a^)'-^a;2fc+2g-..t^^ ^ (3g) 

"'0 

where k = \m\ + p, and p = ii I — \m\ is even, and p = 1 otherwise. Performing the 
integration, we obtain 

^'^^'^'^^" 7r[(2|m| + 2p+l)!!p ^^^ ' 

We emphasize that the power indices / + |m| + p + 3 (in absolute value) in rotating Kerr 
spacetimes are smaller than the corresponding power indices (the well known 2/ + 3) in 
spherically symmetric Schwarzschild spacetimes. (There is an equality only for the / = 
\m\, |m| + l modes). This implies a s/ower decay of perturbations in rotating Kerr spacetimes. 
From Eq. ( ^OD it is easy to see that the time scale tc at which the late-time tail of rotating 
gravitational collapse is considerably different from the corresponding tail of non rotating 
collapse (for the / > |m| + 1 modes) is tc = yy'/a, where y' is roughly the average location 
of the initial pulse. 

We further consider the asymptotic behaviour of the field at future null infinity scri^. 
Using Eq. (|32|) , togather with the representation Eq. (^) for the spheroidal wave functions, 
one finds that the behaviour of the / mode at the asymptotic region of null infinity is 
dominated by the following effective Green's function: 

rC(,,'.f)- f M2"^(2A; + l)!(-l)(^+^)/\ ,^, 
'^' '^".^^ vrA;![(2A; + l)!!]2 ^ 

k=\m\+'p L\ / J 

r—iao 

CkiYii9)Y*{9') / (acj)'-^w'=+^e-^"(*~^)rfw , (41) 

Jo 

where p is defined above. Performing the integration we obtain 

c.. ,._ f M2-n-l)('+^+^)/^(2fc + !)!(/ + !)! , ^, 

' ^ ' ' ^ " M., nkim + iy.ir y 

CMYi{e)Y:{e')a'-''u-^'+'^ . (42) 



14 



Finally, we consider the behaviour of the field at the black-hole outer-horizon Hj^. Fol- 
lowing an analysis similar to the one presented in Sec. [IV C| , we find that the asymptotic 
behaviour of the / mode (where / > |m|) at the black- hole outer horizon if+ is dominated 
by the following effective Green's function: 

^'^'^'"^'^^"^'^ vr[(2|m|+2p+l)!!p ^ 



VI. SUMMARY OF RESULTS 

The asymptotic late-time evolution of a test scalar field on a realistic Kerr background 
is characterized by inverse power-law tails at the three asymptotic regions: timelike infinity 
i+, future null infinity scrij^, and the black-hole outer- horizon if+ (where the power-law 
behaviour is multiplied by an oscillatory term, caused by the dragging of reference frames 
at the event horizon). This late-time behaviour is similar to the the late-time evolution of 
massless fields on a s]9/ienca//^symmetric Schwarzschild background, originally analyzed by 
Price 0. 

However, there is one important feature (besides the oscillatory behaviour along the 
black-hole horizon) which is unique to rotating collapse: Active coupling of different modes. 
Physically, this phenomena is caused by the dragging of reference frames, due to the black- 
hole (or star's) rotation (this phenomena is absent in the non-rotating a = case). Hence, 
the late-time tails of realistic rotating collapse are dominated by a mode, which, in general, 
is different from the initial one. 

As shown in our analysis, there are various cases, which result in different types of 
asymptotic late-time evolutions (different damping exponents). These distinct cases are 
characterized by the values of /* and |r7i|, and the parity of /* — \m\. The dominant modes 
at asymptotic late-times and the values of the corresponding damping exponents (for pure 
initial data, characterized by a particular spherical harmonic function Yjl") are summarized 
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in Tables | (for /* > \m\ + 2) and || (for /* = \m\, \m\ + 1). These analytical results should 



be compared with the numerical results of [|T^ - we find a complete agreement between the 
two, as can be seen in Table |IT1[ 

In summary, the rotation of a black hole (or a star) results in an active coupling of 
different modes. Thus, new modes, which are different from the original one, would be 
generated during a rotating collapse and may play an important role at asymptotic late 
times. Hence, the radiative tail of rotating gravitational collapse decays slower than the 
corresponding tail of a non-rotating collapse. In general, a rotating Kerr black hole becomes 
'bald' slower than a spherically-symmetric Schwarzschild black hole. 
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APPENDIX A: COUPLING OF THE REAL AND IMAGINARY PARTS OF * 

In general, the \l/'^ functions are complex quantities. Owing to the purely real character 
of the Bi{r, 9) coefficient, the real and imaginary parts of each \l/ are decoupled for axially- 
symmetric {m = 0) modes (where B2 vanishes and the Green's function is purely real). 
However, for non-axially symmetric {m 7^ 0) modes, the imaginary coefficient B2{r) mixes 
the real and imaginary parts of the same mode. Hence, the real part of a particular mode will 
always generate the corresponding imaginary part of the same mode, and vice versa (note, 
however, that the coefficient B2{r) is ^—independent. Hence, contrasted with the Bi{r,6) 
coefficient, it cannot mix nearest neighbors modes). We now assume that the initial data 
for \& is purely real, and (without loss of generality) we assume that its angular distribution 
is characterized as before by a particular spherical harmonic function YiT. 

We consider first the asymptotic behaviour of the imaginary part of \1/™' at timelike 
infinity i_|_. Following an analysis similar to the one presented in Sec. |1V A] , one finds that 
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the asymptotic late-time behaviour of the imaginary part of the / mode (where |/* — /| is 
even and I > \m\ > 0) is dominated by the following effective Green's function: 

'^' '^".^S.. -[(2A: + 1)!!P ^''^ 

CuiCui*Yi{d)Y;{e')a'*^'~^H~^'*^'+^^ , (Al) 

where p = if /* — \m\ is even and p = 1 otherwise. Here, L = I* for / > /* modes and 
L = I for I < I* — 2 modes. Hence, we find that the evolution at timelike infinity i+ of the 
imaginary part of ^ is dominated by the lowest allowed mode. 

We further consider the asymptotic behaviour of the imaginary part of ^ at future null 
infinity scrij^. Following an analysis similar to the one presented in Sec. |1V lj| , one finds 
that the asymptotic late-time behaviour of the imaginary part of the / mode is dominated 
by the following effective Green's function: 



Gf{z,z'-t) 

for / > /* modes, and 

G?{z,z'-t) 



7r/*![(2/* + l)!!]2 

Ci,iYi{e)Y;{e')a'-''u-^'+^'^ , (A2) 



ci ^ ^ _ M(-l)('*+'+2)/22-'(2/ + !)!(/* + 1)! ,;+, 



7r/!p + l)!!]2 
Cu*Yi{e)Y;{e')a''-'u-^''+^'^ , (A3) 



for I < I* — 2 modes. Thus, we find that the evolution at null infinity scri^ of the imaginary 
part of \E' is dominated by the / < /* modes (provided that |/ — /*| is even and / > \m\ > 0). 
Finally, we note that the Green's function Eq. (^) is already a complex function for 
non-axially symmetric {m ^ 0) modes. Thus, the real and imaginary parts of the same mode 
are coupled by both the Bi{r,9) and B2{r) coefficients. Since the coupling can occur by 
the ^—dependent coefficient Bi{r,6) (which couples nearest neighbors modes), the asymp- 
totic behaviour at the black-hole outer horizon H^ of the imaginary part of the / mode is 
dominated by the effective Green's function Eq. (pHI). 
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TABLES 

TABLE L Dominant modes and asymptotic damping exponents for I* > \m\ + 2 

asymptotic region dominant mode(s) damping exponent 

timelike infinity |?Ti|+p —{l* + \m\+p+l) 

null infinity \m\+ p < I < I* — 2 —I* 



outer horizon ?7i +p 


-{l* + \m\+p+l) 


TABLE IL Dominant mode and asymptotic damping exponents for /* 


= \m\, \m\ + 1 


asymptotic region dominant mode 


damping exponent 



timelike infinity I* -{21* +3) 

null infinity I* -{I* + 2) 

outer horizon I* -{21* + 3) 

TABLE in. analytical vs. numerical results - the dominant damping exponent at i+ 

initial data analytical result numerical result 

l = m = l -5 -4.87 

l = m = 2 -7 -6.98 

/ = 2,m = -3 ~ -2.9 
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